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^ I Abstract 

We describe the neutron polarimeter NP0L3 for the measurement of polarization 
transfer observables Dij with a typical high resolution of ~300 keV at T n ~ 200 
MeV. The NPOL3 system consists of three planes of neutron detectors. The first 
two planes for neutron polarization analysis are made of 20 sets of one-dimensional 
position-sensitive plastic scintillation counters with a size of 100 cm x 10 cm x 5 cm, 
and they cover the area of 100 x 100 cm 2 . The last plane for detecting doubly scat- 
tered neutrons or recoiled protons is made of the two-dimensional position-sensitive 
liquid scintillation counter with a size of 100 cm x 100 cm x 10 cm. The effective 
analyzing powers A y]e g and double scattering efficiencies eo.s. were measured by us- 
ing the three kinds of polarized neutrons from the 2 H(p, n)pp, 6 Li(p, n) 6 Be(g.s.), and 
12 C(p,n) 12 N(g.s.) reactions at T p = 198 MeV. The performance of NPOL3 defined 
as e D.s.(^4j/;eff) 2 are similar to that of the Indiana Neutron POLarimeter (INPOL) 
by taking into account for the counter configuration difference between these two 
neutron polarimeters. 
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1 Introduction 



The polarization transfer observables for the charge-exchange (p, n) reac- 
tion at intermediate energies T p > 100 MeV provides a potentially rich source 
of information not only on nuclear responses but also on effective interac- 
tions, and their extensive studies have been performed at the Los Alamos 
Meson Physics Facility (LAPMF), the Indiana University Cyclotron Facility 
(IUCF), and the Research Center for Nuclear Physics (RCNP). An example 
of these studies is the search for a pionic enhancement in nuclei via measure- 
ments of a complete set of Dij for (p,n) quasielastic scatterings (QES). The 
spin-longitudinal (pionic) response function Rl is theoretically expected to be 
enhanced relative to the spin-transverse response function 0, 0| • The en- 
hancement of Rl is attributed to the collectivity induced by the attraction of 
the one-pion exchange interaction, and has aroused much interest in connec- 
tion with both the precursor phenomena of t he p ion condensation 
and the pion excess in the nucleus 0, S H S Efl ■ Surprisingly, the experimen- 
tally extracted Rl/ Rt ratios are less than or equal to unity 0, 0], which 
contradicts the theoretical predictions of the enhanced Rl and the quenched 



Rt- Recent analyses of the QES data [13J, |1J| show the pionic enhancement 
in the spin-longitudinal cross section which well represents the Rl- The dis- 
crepancy in Rl/Rt might be due to the effects of the medium modifications 
of the effective NN interaction. These effects could be studied by measuring 
a complete set of for stretched states 0, 0] , which requires a relatively 
better energy resolution of ~500 keV compared with the energy resolution of 
2-3 MeV in the QES measurements. 

Another requirement for high resolution is for the measurement of the Gamow- 
Teller (GT) unit cross section cr GT . Recently Yako et al. 0] applied multipole 
decomposition analysis both to their 90 Zr(n,p) data and to the 90 Zr(p, n) data 
in Refs. 003, and obtained the GT quenching factor Q = 0.88 ± 0.03 ± 
0.16. The first uncertainty contains the uncertainties both of the MDA and of 
the estimation of the isovector spin-monopole (IVSM) contributions, and the 
second uncertainty originates from the uncertainty of <3"gt- This large Q value 
clearly indicates that the configuration mixing is the main mechanism of the 
quenching and thus the admixture of the Ah states into the low-lying states 
plays a minor role. However, a relatively large uncertainty of <jgt makes it 
difficult to draw a definite conclusion from this Q value. In principle, a precise 
<jgt value can be obtained by measuring the cross section at 0° for the ground 
or low-lying discrete GT state whose B(GT) value is known by the /3-decay 
[2oj | . In practice, such a measurement is hampered by a poor energy resolution 
of the neutron time-of-flight (TOF) system. Actually, for the TOF facility at 
RCNP, the energy resolution is limited to be ~1.6 MeV at T p ~ 300 MeV 



due to the thick neutron counter thickness of NPOL2 [2l( compared to the 



relatively short flight pass length 100 m |22j. This poor energy resolution does 
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not allow such direct determination of <7gt- 

The demand for higher resolution necessary for nuclear spectroscopy led to the 
design and construction of the new neutron detector and polarimeter NPOL3. 
The neutron detector should be designed so that the final energy resolution 
is better than ~500 keV in order to resolve GT and stretched states from 
their neighboring states. This can be achieved by using the neutron detector 
material thinner than the NPOL2. Furthermore, the NPOL3 is designed to 
realize the Dij measurement which is important to determine the effective NN 
interaction. 

In Section 2, we present the results of the simulation performed to design the 
NPOL3 system. In Section 3, we will describe the NPOL3 and its performance 
for the neutron detector. Sections 4 and 5 are devoted to the calibration and 
optimization procedures of the NPOL3 for the neutron polarimeter. In Sec- 
tions 6 and 7, we will discuss the results of the calibration. A summary is 
given in Section 8. 



2 Improvements of time and energy resolution 

2.1 Time and energy resolution 

At intermediate energies where neutron kinetic energies T n are determined by 
the TOF technique over a long flight path length of ~100 m, a large volume 
of the scintillator is required in order to achieve a sufficient neutron detection 
efficiency. The energy resolution AT n by the TOF technique is related to the 
uncertainties both of timing At/t and flight path length AL/L as 



with 




where m n is the neutron mass. At is the uncertainty of the flight time t which 
originates from the time spread of the incident proton beam and the time res- 
olution of the neutron counter. AL is the uncertainty of the flight path length 
L which comes from the thickness of the counter. The charged particle scat- 
tered by an incident neutron passes the neutron counter and optical photons 
are generated along the path. The number of photons Ar photon is proportional 
to the thickness AL and the time resolution At becomes better by increasing 
-Nphoton- However, the thickness AL is directly related to the energy resolu- 
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tion AT n as is seen in Eq. (1). Thus the AL dependence on AT n is rather 
complicated since At is also the function of AL. 



2.2 Monte-Carlo simulation with GEANT4 



We performed the Monte-Carlo simulation in order to describe the energy res- 
olution T n quantitatively as a function of AL. The computer library geant4 
[iil ] was used to simulate the neutron-induced reactions in scintillator mate- 
rial, the generation of optical photons, and their propagation. The scintillator 
material is assumed to be the plastic scintillator Saint-Gobain |24[ Bicron 
BC408 with a hydrogen to carbon ratio of H/C = 1.1. The size of the scintil- 
lator is 100 cm x 10 cm x AL cm, and its configuration is the one-dimensional 
position-sensitive detector (hodoscope) with a thickness of AL cm (see Fig. 1). 
The length of 100 cm is same as that of the NPOL2, and the width is fixed 
to be 10 cm. At both ends (L:left and R:right) of the scintillator, the timing 
information of arrived optical photons is accumulated. The neutron flight time 
t can be deduced by using timing information, tz and tn, at both ends as 

where t^F corresponds to the timing information for an accelerator RF signal 
in the experiment. In practice, photo-multiplier tubes (PMTs) are mounted 
at both ends and optical photons are converted to electric signals. Each anode 
signal from PMT is fed to a constant fraction discriminator (CFD) to create 
a fast logic timing signal. We simulated the operation of the CFD. The upper 
panel of Fig. 2 is the pulse-height spectrum generated by our simulation, and 
the lower panel is the simulated CFD timing spectrum. The timing informa- 
tion is determined from the zero-crossing point of the CFD spectrum. Thus 
the resulting t^ and tn are almost independent of an arbitrarily set discrimi- 
nation level of the CFD. The time spread of the zero-crossing points in Fig. 2 
corresponds to the counter-thickness (AL) effect to At. 

Figure 3 represents the expected energy resolutions AT n evaluated from the 
simulation results as a function of AL. The filled circles and the filled boxes 
correspond to the results for the measurement at T n = 200 and 300 MeV, 
respectively. In the simulation, the flight path length L is 100 m and the time 
spread of the proton beam is assumed to be 200 ps FWHM which is a typical 
value of the beam from the RCNP Ring cyclotron. The AT n values clearly 
depend on AL for both T n = 200 and 300 MeV. It is found that the counter 
thickness AL should be thinner than 5 cm in order to achieve the required 
final resolution of ~500 keV at T n = 300 MeV. Thus we have employed the 
hodoscope counters with a thickness of 5cm in the NPOL3. 
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2.3 Comparison with experimental results 

We have constructed the one- dimensional plastic scintillation counters with 
a size of 100 cm x 10 cm x 5cm. With these counters, we measured the 
12 C(p, n) reaction at 9\ a b = 0° and incident beam energies of T p = 198 and 
295 MeV. The target thickness is 38 mg/cm 2 whose contribution to the final 
energy resolution is negligible small. The results are shown in Fig. 5. The 
energy resolutions were evaluated by fitting the peak of the 12 N state to the 
standard hyper-Gaussian and they are about 280 and 500 keV FWHM for T p 
= 198 and 295 MeV data, respectively. These energy resolutions are consistent 
with the simulation results shown in Fig. 3, and they satisfy our requirements. 
These counters have been used in the NPOL3 which will discussed in the next 
section. 



3 Neutron detector NPOL3 and its performance 

3. 1 Neutron detector NPOL3 

The neutron detector NPOL3 has two main neutron detector planes, HD1 
and HD2. Each plane consists of 10 stacked 100 cm x 10 cm x 5cm plastic 
scintillator bars made of Bicron BC408. Two thin plastic scintillator planes 
CPV and CPD, used to identify charged particles, complete the NPOL3 which 
is schematically shown in Fig. 4. Additional two-dimensional neutron counter 
NC is used in polarimetry mode, which will be discussed in Sec. 4. A scintillator 
in HD1 or HD2 is viewed by the fast Hamamatsu 25] H2431 (HD1) or H1949 
(HD2) PMTs attached at each end. The one-dimensional position is deduced 
from the fast timing information derived from PMT. The CPV and CPD are 
made of three sets of plastic scintillators (BC408) with a size of 102 cm x 35 
cm x 0.5 cm. The scintillator is viewed from both sides (U:up and D:down) by 
the Hamamatsu H7195 PMT through fish-tail shape light guides. The CPV is 
used to veto charged particles entered into the NPOL3, and the CPD is used 
to detect the recoiled charged particles. 

In the neutron detector mode of NPOL3, a neutron is detected by either 
neutron detector plane HD1 or HD2. Furthermore we have required that the 
recoiled charged particle should be detected by the following detector plane. 
This means that we measure neutrons in the polarimetry (double-scattering) 
mode of the NPOL3 system. This procedure has been applied in the data 
analysis of the INPOL system j^], and it is useful to improve the final energy 
resolution as well as to reduce the low-energy tail component. Thus we have 
applied this procedure in our data analysis. 
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3.2 Light output calibration 



The neutron detection efficiency depends on the threshold level £ t h of the 
scintillator light output. The total light output £ is constructed by taking the 
geometric mean of the PMT outputs at both ends, Cl and jCr, as 



The light output C is calibrated by using 4.4 MeV 7-rays from a 241 Am- 9 Be 
neutron source and cosmic rays (mostly fi^). In the cosmic-ray measurement, 
10-bar hits in either HD1 or HD2 were required. Thus horizontal and zenith 
angles can be measured in HD1 and HD2, respectively, on the assumption 
that cosmic rays pass in straight lines through the detector plane. The £ for 
cosmic rays has been corrected for these angles in order to evaluate the C in 
case cosmic rays pass the scintillator of 10 cm in thickness. 

The left and right panels of Fig. 6 show the C spectra for 4.4 MeV 7-rays and 
cosmic rays, respectively. We can clearly see the Compton edge (4.2 MeV cc ) for 
the 7-ray spectrum and the sharp peak with a Landau tail for the cosmic-ray 
spectrum. The histograms represent the results of the Monte-Carlo simula- 
tions with geant4 described in Sec. 2. The simulations successfully repro- 
duce the main components of the Compton edge and the sharp peak in 7-ray 
and cosmic-ray spectra, respectively. Furthermore they well describe the one- 
photon escape peak at ~3.5 MeV ee in the 7-ray spectrum and the Landau 
tail in the cosmic-ray spectrum. Note that the discrepancy between the data 
and the simulation result for 4.4 MeV 7-rays in the low C region is due to the 
contributions from 7-rays other than 4.4 MeV 7-rays from the neutron source 
which are not taken into account in the simulations. 

3.3 Neutron detection efficiency 

The laboratory differential cross section da^/dVL is related to the number of 
measured neutrons N„ as 



where N p is the number of incident protons, p is the target density, AQ is the 
solid angle, e is the intrinsic neutron detection efficiency, T is the transmission 
factor along the neutron flight path, and fu ve is the experimental live time. 

We measured the the product eT by measuring the neutron yield from the 0° 
7 Li(p, n) 7 Be(g.s. + 0.43 MeV) reaction which has a constant center of mass 




(4) 



daiab N n 



(5) 



dVL " N pP AneTf hvc ' 
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(cm.) cross section of a r rn (0°) = 27.0 ± 0.8 mb/sr at the bombarding energy 
range of T p = 80-795 MeV j27J. The measurements were performed at T p = 198 
and 295 MeV by using the enriched 7 Li target with a thickness of 54 mg/cm 2 . 
The integrated beam current N p was measured by using an electrically isolated 
graphite beam stop (Faraday cup) connected to a current integrator. The live 
time /ii ve was typically 90%. 

The finally obtained neutron detection efficiencies (eT) are shown in Fig. 7 
as a function of the light output threshold Ah- The systematic uncertainty 
is estimated to be about 6% by taking into account of the uncertainties both 
of the 7 Li(p, n) cross section and of the 7 Li target thickness. It is found that 
the detection efficiencies are almost independent of the neutron kinetic energy 
with a value of ~0.017 at £ t h = 5 MeV ee - 



4 Neutron polarimeter NPOL3 and calibration procedure 

4-1 Neutron polarimeter NPOL3 

In the polarimetry mode of NPOL3, the neutron detector planes, HD1 and 
HD2, serve as neutron polarization analyzers, and the following two-dimensional 
position-sensitive neutron counter NC acts as a catcher of doubly scattered 
neutrons or recoiled protons. Neutron polarization is determined from the 
asymmetry of the n + p events whose analyzing powers can be rigorously ob- 
tained by using the phase-shift analysis of the NN scattering. The n+p events 
are kinematically resolved from the other events by using time, position, and 
pulse-height information from both analyzer and catcher planes. Note that 
this kinematical selection significantly reduces the background events from 
the wraparound of slow neutrons from preceding beam pulses, cosmic rays, 
and the target gamma rays. However, the quasielastic n + C reaction has an 
influence on the determination of the neutron polarization because its kine- 
matical condition is similar to that of the n + p scattering. Thus we have 
measured the effective analyzing powers which include the contributions from 
both n + p and n + C by using the polarized neutron beams. 

4 ■ 2 Effective analyzing powers 

There are two channels, (n,n) and (n,p), in the polarimetry mode. Doubly 
scattered neutrons and recoiled protons are measured by the catcher plane 
NC in (n, n) and (n, p) channels, respectively. The effective analyzing powers, 
A™ s and Ay^ eS , of (n, n) and (n, p) channels can be determined by using the 
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neutron beam with a known polarization p' N . Ay™ s and Ay? eS are deduced by 

^nn ^np 

A y?es = A y;eS = ~r> ( 6 ) 

Pn Pn 

where e nn and e np are the asymmetries measured for (n, n) and (n, p) channels, 
respectively. 



4-3 Polarized neutron beams 



The first source of polarized neutrons is the 0° 2 H(p, n)pp reaction at T p = 
198 MeV. The Dll value at T p = 198 MeV was obtained by using the data 
in Ref. j^. Because this reaction is mainly a Gamow- Teller transition, the 
polarization transfer coefficients D{j satisfy ^| 

The neutron polarization p' N is given by using Dnn(0°) as 

p' N = p N D NN (0°), (8) 
where p^ is the polarization of the incident proton beam. 

The second source is the 0° 6 Li(p, ra) 6 Be(g.s.) reaction at T p = 198 MeV. This 
reaction is also a GT transition, and its D NN value was measured at T p = 
200 MeV by Taddeucci 29]. The result can be used as a double-check of the 
calibration performed by using the neutrons from 2 H(p, n)pp. 

The third source is the 0° 12 C(p, n) 12 N(g.s.) reaction at T p = 198 MeV. The 
Dnn{0°) value of this GT transition was also measured at T p = 200 MeV by 
Taddeucci j2^|. Unfortunately, its uncertainty is relatively large compared with 
that of the 2 H(p, n) or 6 Li(jo, n) 6 Be(g.s.) reaction. Thus the final accuracy of 
Ay- e s is limited by the uncertainty of the neutron beam polarization originating 
from the uncertainty of Dnn{0°), Therefore we have used the following method 
which does not require to know the Da(0°) values. 

The GT transition is the spin-flip and unnatural-parity transition, and its po- 
larization transfer coefficients satisfy the relation in Eq. (7). Here we assume 
that two kinds of polarized proton beams are available; one has pure longi- 
tudinal polarization pi and the other has pure normal polarization p^- The 
neutron polarizations at 0° become p' L = PlDll(0°) and p' N = pnDnn(0°) for 
the beam polarizations of pl and p^, respectively. Then the resulting asym- 
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metries measured by a neutron polarimeter are 



e L = p' L A y;cS = p L D LL (0 o )A y . cS , 

e N = PN^-y^S = PND NN (0°)A yie s. 



(9) 



By using Eqs. (7) and (9), Ay- e s can be expressed as 



\Pl PnJ 



(10) 



Thus the A y . c g value can be obtained without knowing a priori the values of 
Du(0°), and we have applied this technique to the 12 C(p, n) 12 N(g.s.) reaction. 

4-4 Experimental conditions 

We used a deuterated polyethylene (CD 2 ) target with a thickness of 228 
mg/cm 2 as a deuteron target. Neutrons from the 2 H(p, n)pp reaction were 
clearly separated from those from the 12 C(p, nX) reaction since the former 
reaction Q value of —2.2 MeV is significantly smaller than the latter Q value 
of > —18.1 MeV. We also used 99% enriched 6 Li and natural C targets with 
thicknesses of 181 and 38 mg/cm 2 , respectively. The proton beam intensity 
and its polarization are typically 80 nA and 0.70, respectively. 

4-5 Sector methods 

We have adopted the sector method in order to obtain the asymmetry. In 
this method, the plane (9, <f>) defined by the scattering polar angle 9 and the 
azimuthal angle is divided into four sectors: Left, Right, Up, and Down, as 
is shown in Fig. 8. Each sector is defined by scattering angle limits 9 mm and 
Qmax an( j j- ne az j mu thal half-angle $. 

The asymmetry e is calculated by using the event numbers in these sectors. 
For example, the relation between the left-right asymmetry e'^ and the event 
numbers is as follows. The azimuthal distribution of scattered particles can be 
described as 



where we suppress the 9 dependence for simplicity. Then the event numbers 
scattered to left and right become 



(11) 




N 



7T 




(12) 




N 




7T 
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where N is the event number accepted in the 9 bin, and the acceptance is 
assumed to be constant. Thus the asymmetry = p' N A y . ef[ is found to be 



e> - Yl ~ Yr $ (13) 

% -y L + ^sin$- {l6) 

Note that the sideway polarization p' s is deduced from the up-down asymmetry 
e' s which is obtained by substituting L — > D and R — ► U in Eq. (13). 



5 Optimization of performance of NPOL3 



5.1 Optimization criteria 



The performance of a polarimeter can be measured by its figure-of-merit 
(FOM) value defined as 

FOM = N(A y , c{i ) 2 , (14) 

where N is the total event number Yl + Yr or Yjj + Yd. We can minimize the 
uncertainty of the final result (neutron polarization and Dij) by maximizing 
this FOM value. 

We have applied three types (Types-I, II, and III) of software cuts to the mea- 
sured 9, 0, inter-plane velocity, and pulse height information of the NPOL3. 
The FOM value is maximized in Type-I, whereas the A y . cS is maximized in 
Type-II with keeping the FOM value larger than 40% of the optimum value. 
In Type-Ill, both FOM and A y . e s values have been optimized with keeping 
the FOM value larger than 80% of the optimum value. The Type-I is the best 
choice from the statistical point of view in order to maximize the performance 
of the NPOL3. However, in general, a higher A y . eS value in Type-II or III is 
helpful to reduce the systematic uncertainty of the final result. Thus we have 
obtained these three sets and we will select the optimum set depending on the 
experimental condition. 

In the following, we show the results by using the data of the 2 H(j9, n)pp 
reaction because the total event number is significantly larger than that of 
other two polarized-neutron production reactions. The A y . c g value is also ob- 
tained from the 6 Li(p, n) 6 Be(g.s.) data by applying the same software cuts. 
The mean neutron energy is about 193 MeV for both the the 2 H(p, n)pp 
and 6 Li(p, n) 6 Be(g.s.) reactions. Thus we can use the results for a double 
check of the calibration. The same software cuts are also applied to the 
12 C(j9, n) 12 N(g.s.) data. The mean neutron energy of ~180 MeV is signifi- 
cantly low, therefore, the result can be used to discuss the neutron energy 
dependence of A y - C s. 
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5.2 Optimization o/$ 



Figure 9 shows the A y . e g and FOM values in Type-I as a function of the 
azimuthal half- angle $. The FOM values are normalized to a maximum FOM 
value of 1. The A y;e s values are insensitive to the choice of $, whereas the 
FOM values take a maximum at $ ~ 67°. This angle is consistent with the 



simple estimation of the optimum $ = 66.8° |2jj. Because the $ dependence 



is almost independent of the software cut types, we set $ = 66.8° in all cases. 



5.3 Optimization of 9 window 



The event selection on 9 is necessary to prevent dilution of A y . e g by events 
with small or negative analyzing powers. Figure 10 shows the scattered-particle 
distribution as a function of 9 in Type-I. The results in other Types-II and 
III are also identical. The double scattering events with a scattering angle up 
to ~ 60° are accepted by the NPOL3. The analyzing power A y {9) of the free 
n+p scattering is positive only up to 45° and 41° for (n, n) and (n,p) channels, 
respectively, at T n = 193 MeV. Thus part of the events should be rejected on 
the basis of small A y as described below. 

Figure 11 shows the contour plots of the FOM values in Type-I as functions 
of both 9 mm and # max . Figure 12 represents the A y . e s and FOM values as a 
function of either 9 mm or # max . In both figures, the FOM values are normalized 
to a maximum FOM value of 1. The FOM values strongly depend on both 9 mm 
and # max ; whereas the A y . e g values weakly depend on them except for # max in 
the (n, n) channel. In cases of Types-II and III, the optimal (9 mm , # max ) values 
are selected on basis of the criteria described in Sec. 5.1. In Table 1, we present 
the finally chosen (9 min , # max ) values. 



5.4 Optimization of R v window 



The kinematical selection of the n + p events is performed by using the kine- 
matical quantity R v defined as 

Rv = — , (15) 

VNN 

where v is the measured particle velocity between the analyzer and catcher 
planes and vjyjy is the predicted velocity based on the NN kinematics. The 
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vnn value is given by 



2cos6JT n /(T n + 2m N ) 

V NN = C- r, 16 

1 - cos 9 T n j [T n + m N ) 

where c is the light velocity and is the nucleon mass. Figure 13 shows 
the scattered-particle distribution in Type-I as a function of R v . The peaks at 
R v ~ 1 correspond to the events from the p(n,n)p or p(n,p)n reaction. The 
shoulders at < 1 are due to the contributions from quasielastic scattering, 
inelastic scattering, and wraparounds. 

Figure 14 shows the contour plots of the FOM values as functions of both 
i?™ in and i?™ ax . Figure 15 represents the A y;e g and FOM values as a function 
of either i2™ m or i?™ ax . The FOM values are normalized to a maximum FOM 
value of 1 in both figures. It is found that the FOM values are almost insensitive 
to i?™ ax for i?™ ax > 1.3. However, i?™ ax has been limited to be less than 1.3 
to eliminated the contribution from 7-ray events with R v > 1.5. In Types-II 
and III, the optimal (-R™ m , -R™ ax ) values are selected on basis of the criteria 
described in Sec. 5.1. The final (-R™ in , -R™ ax ) values chosen are listed in Table 1. 



5.5 Optimization o/£ th 



Figure 16 shows the A y . eS and FOM values in Type-I as a function of C t h- 
The FOM values monotonously decrease as increasing C th , whereas the A y . e s 
values are almost independent of £ t h- In all software cut types, we have chosen 
£ t h = 4 MeV ee to ensure reproducibility in the offline analysis rather than just 
accept the hardware threshold of ~2 MeV ee . 



6 Performance of polarimeter NPOL3 



6. 1 Effective analyzing powers 



The effective analyzing powers A y . ea obtained from the 2 H(p, n)pp data by 
applying previously described software cuts are summarized in Table 2. The 
A y]e s values deduced from the 6 Li(p, n) 6 Be(g.s.) data are also listed in Table 2. 
The results of the 2 H(p, n)pp and 6 Li(p, n) 6 Be(g.s.) data are consistent with 
each other within their statistical uncertainties, which suggests the reliability 
of the present calibration procedure. 

The A y . ef [ values deduced from the 12 C(p, n) 12 N(g.s.) data are plotted in 
Fig. 17, and they are also listed in Table 2. The results of the 2 H(p, n)pp 
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data are also displayed in Fig. 17. The A y - e g values have been estimated as a 
function of T n on the basis of the free n + p analyzing powers derived from 
the NN phase shift analysis j^. In the calculations, the (9,4>) cuts and the 
counter geometries have been properly taken into account. The solid curves 
in Fig. 17 show the results of the calculations. The measured A y;e s values are 
different from the calculated values by factors of 0.73-0.79 and 1.16-1.29 for 
(ft, n) and (n,p) channels, respectively, depending on the event cut types. The 
difference of A y . e g is mainly due to the contributions from the quasielastic 
n + C events. The dashed curves represent the calculated A y . eS values normal- 
ized to reproduced the experimental results. These normalized values will be 
used in future data analysis in case there is no appropriate calibration result. 



6.2 Double scattering efficiency 

The double scattering efficiencies (cd.s.) of the NPOL3 were deduced from 
the data of the 2 H(p, n)pp and 12 C(p, n) 12 N(g.s.) reactions. The results are 
listed in Table 3. The 6d.s. values in both (ft, n) and (n,p) channels are almost 
constant in the present neutron energy region. 



7 Discussion 



7.1 9 dependence of A y . eS 

In the sector method, the events within (9 mm , # max ) are integrated in order to 
evaluate the A y . C Q values in the calibration. Thus the results can be considered 
as the mean effective analyzing powers in the range of 9 = # mm -# max . We 
deduced the A y . eS values as a function of 9 in order to check whether the 
n + p events are properly selected in our analysis. As a typical example, the 
results in Type-Ill are shown in Fig. 18. The solid curves display the angular 
distributions of A y (9) of the free n + p scattering. The measured angular 
distributions are very similar to those of the free n+p scattering. However the 
magnitudes are different mainly due to the contributions from the quasielastic 
n + C events. The dashed curves are the results normalized to reproduce 
the measured values, and the normalization factors are 0.7 and 1.2 for (n, n) 
and (n, p) channels, respectively. The agreement of the angular distributions 
supports the proper kinematical selection of the n + p events in our analysis. 
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7.2 Comparison with INPOL system 



Here we compare the performance of the NPOL3 system with that of the 
INPOL system [26| which has been used in the same energy region. For the 
(n, n) channel, the NPOL3 has ~l/2 volume for both analyzer and catcher 
planes compared with the INPOL. Thus the double scattering efficiency cd.s 
and the corresponding FOM value of the NPOL3 is expected to be ~l/4 of 
those of the INPOL. The FOM value of the NPOL3 is 0.28 x 1(T 4 in Type- 
I at T n = 193 MeV, which is about 1/4 of that of the INPOL reported as 
0.92 x 10~ 4 at T n = 194 MeV. For the (n,p) channel, the FOM value of the 
NPOL3 is expected to be similar to that of the INPOL because their volume 
of the analyzer planes is similar with each other. In fact, the FOM value of 
0.89 x 10~ 4 in Type-I is close to that of 0.81 x 10~ 4 of the INPOL. Thus we 
conclude that the performance of the NPOL3 is very similar to that of the 
INPOL by taking into account of the difference of the thicknesses of both 
analyzer and catcher planes. 



8 Summary 



The high resolution neutron polarimeter NPOL3 has been constructed and 
developed to measure polarization transfer observables Dij for (p, n) reactions 
at intermediate energies around T p = 200 MeV. The NPOL3 system can mea- 
sure the normal as well as sideways components of the neutron polarization 
simultaneously. The other longitudinal component can be measured by using 
the neutron spin rotation magnet described in Ref. ji^. Thus we can perform 
the measurement of a complete set of with a high resolution of ~300 keV 
by using the NPOL3. 

The performance of the NPOL3 system was studied by using the polarized 
neutrons from 2 H(p, n)pp, 6 Li(p, n) 6 Be(g.s.), and l2 C(p, n) 12 N(g.s.) reactions 
at T p = 198 MeV. The effective analyzing powers in Type-I are 0.33 and 0.12 
for (ft, n) and (n,p) channels, respectively. The performance is comparable to 
that of the INPOL system at the same energy. 
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Neutron 



Fig. 1. Illustration of the Monte-Carlo simulation with geant4. The neutron beam 
with T n = 200 or 300 MeV bombard a plastic scintillator with a size of 100 cm x 10 
cm x AL cm. The optical photons generated by recoiled charge particles propagate 
in the scintillator and are detected by PMTs at both ends. The time information, 
tL and tR, deduced from PMT signals are used to evaluate the time of flight t and 
its width At. 




Time (ns) 

Fig. 2. The pulse-height (upper panel) and CFD timing (lower panel) spectra 
generated by the Monte-Carlo simulation with geant4. See text for details. 
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Fig. 3. The simulation results for the energy resolution AT n as a function of the 
counter thickness AL. The filled circles and the filled boxes represent the results for 
the measurement at T„ = 200 and 300 MeV, respectively, as explained in the text. 




Fig. 4. Schematic (left) and top (right) views of the NPOL3 system. In the po- 
larimetry mode of NPOL3, HD1 and HD2 are the analyzer planes while NC is the 
catcher plane. Thin plastic scintillator planes are used to veto (CPV) or identify 
(CPD) charged particles. 
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Excitation energy of "N (MeV) 



Fig. 5. Measurements of the 12 C(p, n) reaction at T p = 198 MeV (upper panel) and 
297 (lower panel) MeV with 5 cm-thick plastic scintillation counters. The curve of 
each panel shows the reproduction of the 12 N(g.s.) state with the hyper-Gaussian 
peak. 




1 2 3 4 5 6 15 20 



Light output £ (MeV ee ) 

Fig. 6. Light output spectra for 4.4 MeV 7-rays from a 241 Am- 9 Be neutron source 
(left panel) and cosmic rays (right panel). The histograms represent the results of 
the Monte-Carlo simulations with geant4. See text for details. 
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Fig. 7. Neutron detection efficiencies of 5-cm thick plastic scintillation counters as a 
function of detection threshold C t h- The filled circles and the filled boxes represent 
the results at T n = 190 and 291 MeV, respectively. 
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Fig. 9. Effective analyzing powers A y]e s and FOM as a function of azimuthal 
half-angle <£. The upper and lower panels correspond to the results for (n, n) and 
(n,p) channels, respectively. 




9 (deg) 

Fig. 10. Distributions of the doubly scattered particles as a function of scattering 
polar angle 9. The left and right panels correspond to the results for (n, n) and 
(n,p) channels, respectively. 
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5 10 15 20 25 5 10 15 20 25 30 

min (deg) 

Fig. 11. FOM of NPOL3 as functions of lower and upper limits of the scattering 
polar angle window (0 mm , max ) . The left and right panels correspond to the re- 
sults for (n, n) and (n, p) channels, respectively. The selected limits in Type-I are 
indicated by the filled circles. 



| l l l l | l l l l | l l l l | — l l l l 
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6 min (deg) max (deg) 

Fig. 12. Effective analyzing powers A y . e g and FOM as a function of either lower (left 
panels) or upper (right panels) limit of the scattering polar angle 6. The upper and 
lower panels correspond to the results for (n, n) and (n,p) channels, respectively. 
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Fig. 13. Distributions of the doubly scattered particles as a function of velocity 
ratio R v . The left and right panels correspond to the results for (n, n) and (n, p) 
channels, respectively. 



S 1.4 h 

B » 




Fig. 14. FOM of NPOL3 as functions of lower and upper limits of the velocity ratio 
window (i?" im , i?™ ax ). The left and right panels correspond to the results for (n,n) 
and (n,p) channels, respectively. The selected limits in Type-I are indicated by the 
filled circles. 
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Fig. 15. Effective analyzing powers A y . e f[ and FOM as a function of either lower 
(left panels) or upper (right panels) limit of the velocity ratio R v . The upper and 
lower panels correspond to the results for (ft,n) and (n,p) channels, respectively. 
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Fig. 16. Effective analyzing powers ^eff and FOM as a function of detection 
threshold C t h- The upper and lower panels correspond to the results for (ft, n) and 
(n,p) channels, respectively. 
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Fig. 17. Effective analyzing powers of the NPOL3 system for (n, n) (filled circles) 
and (n,p) (filled boxes) channels, respectively, as a function of neutron energy. The 
top, middle, and bottom panels correspond to the results in the software cut Types-I, 
II, and III, respectively. The solid curves represent the calculations based on the 
free n + p observables. The dashed curves represent the calculations normalized to 
reproduce the experimental data. 
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6 (deg) 



Fig. 18. Angular distributions of effective analyzing powers of the NPOL3 system 
for (n, n) (filled circles) and (n,p) (filled boxes) channels, respectively. The solid 
curves represent the calculations based on the free n + p observables whereas the 
dashed curve represent the normalized results to reproduce the experimental data. 
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Type-I Type-II Type-Ill 

Quantity FOM optimum A y . e ff optimum Intermediate 





(n, n) 


(n,p) 


(n, n) 


(n,p) 


(re, n) 


(re» 




67° 


67° 


67° 


67° 


67° 


67° 


8 


(8°, 30°) 


(8°, 53°) 


(12°, 24°) 


(20°, 53°) 


(10°, 27°) 


(12°, 53°) 


R v 


(0.90,1.30) 


(0.70, 1.30) 


(0.96,1.30) 


(0.88,1.30) 


(0.92,1.30) 


(0.70,1.30) 


Ah 


4MeV ee 


4MeV ee 


4MeV ee 


4MeV ee 


4MeV ee 


4MeV ee 



Table 1 



Optimized cuts used in the NPOL3 system. 



Type-I Type-II Type-Ill 

Energy (Reaction) FOM optimum A y . e fi optimum Intermediate 





(re, n) 


(n,p) 


(n, n) 


(n,p) 


(re, n) 


(n,p) 


193 MeV ( 2 H(p, n)pp) 


0.333 


0.123 


0.398 


0.154 


0.373 


0.128 




±0.008 


±0.002 


±0.013 


±0.003 


±0.009 


±0.002 


193 MeV ( 6 Li(p,n) 6 Be(g.s.)) 


0.305 


0.123 


0.369 


0.148 


0.341 


0.136 




±0.025 


±0.009 


±0.032 


±0.011 


±0.028 


±0.010 


180 MeV ( 12 C(p,ra) 12 N(g.s.)) 


0.358 


0.105 


0.427 


0.129 


0.398 


0.112 




±0.021 


±0.004 


±0.034 


±0.008 


±0.024 


±0.005 



Table 2 



Effective analyzing powers of the NPOL3 system. 



Type-I Type-II Type-Ill 

Energy FOM optimum ^4y ;c ff optimum Intermediate 

(n,n) (n,p) (n,n) (n,p) (ft,n) (n,p) 

193 MeV 2.5 x 10~ 4 5.9 x 10~ 3 0.9 x 10~ 4 1.7 x 10~ 3 1.8 x 10~ 4 4.9 x 10~ 3 

180 MeV 2.8 x 10~ 4 6.5 x 10~ 3 1.0 x 10~ 4 1.8 x 10" 3 2.0 x 10~ 4 5.4 x 10" 3 
Table 3 

Double scattering efficiencies of the NPOL3 system. 
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We describe the neutron polarimeter NP0L3 for the measurement of polarization 
transfer observables Dij with a typical high resolution of ~300 keV at T n ~ 200 
MeV. The NPOL3 system consists of three planes of neutron detectors. The first 
two planes for neutron polarization analysis are made of 20 sets of one-dimensional 
position-sensitive plastic scintillation counters with a size of 100 cm x 10 cm x 5 cm, 
and they cover the area of 100 x 100 cm 2 . The last plane for detecting doubly scat- 
tered neutrons or recoiled protons is made of the two-dimensional position-sensitive 
liquid scintillation counter with a size of 100 cm x 100 cm x 10 cm. The effective 
analyzing powers A y]e g and double scattering efficiencies eo.s. were measured by us- 
ing the three kinds of polarized neutrons from the 2 H(p, n)pp, 6 Li(p, n) 6 Be(g.s.), and 
12 C(p,n) 12 N(g.s.) reactions at T p = 198 MeV. The performance of NPOL3 defined 
as e D.s.(^4j/;eff) 2 are similar to that of the Indiana Neutron POLarimeter (INPOL) 
by taking into account for the counter configuration difference between these two 
neutron polarimeters. 
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1 Introduction 



The polarization transfer observables for the charge-exchange (p, n) reac- 
tion at intermediate energies T p > 100 MeV provides a potentially rich source 
of information not only on nuclear responses but also on effective interac- 
tions, and their extensive studies have been performed at the Los Alamos 
Meson Physics Facility (LAPMF), the Indiana University Cyclotron Facility 
(IUCF), and the Research Center for Nuclear Physics (RCNP). An example 
of these studies is the search for a pionic enhancement in nuclei via measure- 
ments of a complete set of Dij for (p,n) quasielastic scatterings (QES). The 
spin-longitudinal (pionic) response function Rl is theoretically expected to be 
enhanced relative to the spin-transverse response function 0, 0| • The en- 
hancement of Rl is attributed to the collectivity induced by the attraction of 
the one-pion exchange interaction, and has aroused much interest in connec- 
tion with both the precursor phenomena of t he p ion condensation 
and the pion excess in the nucleus 0, S H S Efl ■ Surprisingly, the experimen- 
tally extracted Rl/ Rt ratios are less than or equal to unity 0, 0], which 
contradicts the theoretical predictions of the enhanced Rl and the quenched 



Rt- Recent analyses of the QES data [13J, |1J| show the pionic enhancement 
in the spin-longitudinal cross section which well represents the Rl- The dis- 
crepancy in Rl/Rt might be due to the effects of the medium modifications 
of the effective NN interaction. These effects could be studied by measuring 
a complete set of for stretched states 0, 0] , which requires a relatively 
better energy resolution of ~500 keV compared with the energy resolution of 
2-3 MeV in the QES measurements. 

Another requirement for high resolution is for the measurement of the Gamow- 
Teller (GT) unit cross section cr GT . Recently Yako et al. 0] applied multipole 
decomposition analysis both to their 90 Zr(n,p) data and to the 90 Zr(p, n) data 
in Refs. 003, and obtained the GT quenching factor Q = 0.88 ± 0.03 ± 
0.16. The first uncertainty contains the uncertainties both of the MDA and of 
the estimation of the isovector spin-monopole (IVSM) contributions, and the 
second uncertainty originates from the uncertainty of <3"gt- This large Q value 
clearly indicates that the configuration mixing is the main mechanism of the 
quenching and thus the admixture of the Ah states into the low-lying states 
plays a minor role. However, a relatively large uncertainty of <jgt makes it 
difficult to draw a definite conclusion from this Q value. In principle, a precise 
<jgt value can be obtained by measuring the cross section at 0° for the ground 
or low-lying discrete GT state whose B(GT) value is known by the /3-decay 
[2oj | . In practice, such a measurement is hampered by a poor energy resolution 
of the neutron time-of-flight (TOF) system. Actually, for the TOF facility at 
RCNP, the energy resolution is limited to be ~1.6 MeV at T p ~ 300 MeV 



due to the thick neutron counter thickness of NPOL2 [2l( compared to the 



relatively short flight pass length 100 m |22j. This poor energy resolution does 
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not allow such direct determination of <7gt- 

The demand for higher resolution necessary for nuclear spectroscopy led to the 
design and construction of the new neutron detector and polarimeter NPOL3. 
The neutron detector should be designed so that the final energy resolution 
is better than ~500 keV in order to resolve GT and stretched states from 
their neighboring states. This can be achieved by using the neutron detector 
material thinner than the NPOL2. Furthermore, the NPOL3 is designed to 
realize the Dij measurement which is important to determine the effective NN 
interaction. 

In Section 2, we present the results of the simulation performed to design the 
NPOL3 system. In Section 3, we will describe the NPOL3 and its performance 
for the neutron detector. Sections 4 and 5 are devoted to the calibration and 
optimization procedures of the NPOL3 for the neutron polarimeter. In Sec- 
tions 6 and 7, we will discuss the results of the calibration. A summary is 
given in Section 8. 



2 Improvements of time and energy resolution 

2.1 Time and energy resolution 

At intermediate energies where neutron kinetic energies T n are determined by 
the TOF technique over a long flight path length of ~100 m, a large volume 
of the scintillator is required in order to achieve a sufficient neutron detection 
efficiency. The energy resolution AT n by the TOF technique is related to the 
uncertainties both of timing At/t and flight path length AL/L as 



with 




where m n is the neutron mass. At is the uncertainty of the flight time t which 
originates from the time spread of the incident proton beam and the time res- 
olution of the neutron counter. AL is the uncertainty of the flight path length 
L which comes from the thickness of the counter. The charged particle scat- 
tered by an incident neutron passes the neutron counter and optical photons 
are generated along the path. The number of photons Ar photon is proportional 
to the thickness AL and the time resolution At becomes better by increasing 
-Nphoton- However, the thickness AL is directly related to the energy resolu- 
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tion AT n as is seen in Eq. (1). Thus the AL dependence on AT n is rather 
complicated since At is also the function of AL. 



2.2 Monte-Carlo simulation with GEANT4 



We performed the Monte-Carlo simulation in order to describe the energy res- 
olution T n quantitatively as a function of AL. The computer library geant4 
[iil ] was used to simulate the neutron-induced reactions in scintillator mate- 
rial, the generation of optical photons, and their propagation. The scintillator 
material is assumed to be the plastic scintillator Saint-Gobain |24[ Bicron 
BC408 with a hydrogen to carbon ratio of H/C = 1.1. The size of the scintil- 
lator is 100 cm x 10 cm x AL cm, and its configuration is the one-dimensional 
position-sensitive detector (hodoscope) with a thickness of AL cm (see Fig. 1). 
The length of 100 cm is same as that of the NPOL2, and the width is fixed 
to be 10 cm. At both ends (L:left and R:right) of the scintillator, the timing 
information of arrived optical photons is accumulated. The neutron flight time 
t can be deduced by using timing information, tz and tn, at both ends as 

where t^F corresponds to the timing information for an accelerator RF signal 
in the experiment. In practice, photo-multiplier tubes (PMTs) are mounted 
at both ends and optical photons are converted to electric signals. Each anode 
signal from PMT is fed to a constant fraction discriminator (CFD) to create 
a fast logic timing signal. We simulated the operation of the CFD. The upper 
panel of Fig. 2 is the pulse-height spectrum generated by our simulation, and 
the lower panel is the simulated CFD timing spectrum. The timing informa- 
tion is determined from the zero-crossing point of the CFD spectrum. Thus 
the resulting t^ and tn are almost independent of an arbitrarily set discrimi- 
nation level of the CFD. The time spread of the zero-crossing points in Fig. 2 
corresponds to the counter-thickness (AL) effect to At. 

Figure 3 represents the expected energy resolutions AT n evaluated from the 
simulation results as a function of AL. The filled circles and the filled boxes 
correspond to the results for the measurement at T n = 200 and 300 MeV, 
respectively. In the simulation, the flight path length L is 100 m and the time 
spread of the proton beam is assumed to be 200 ps FWHM which is a typical 
value of the beam from the RCNP Ring cyclotron. The AT n values clearly 
depend on AL for both T n = 200 and 300 MeV. It is found that the counter 
thickness AL should be thinner than 5 cm in order to achieve the required 
final resolution of ~500 keV at T n = 300 MeV. Thus we have employed the 
hodoscope counters with a thickness of 5cm in the NPOL3. 
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2.3 Comparison with experimental results 

We have constructed the one- dimensional plastic scintillation counters with 
a size of 100 cm x 10 cm x 5cm. With these counters, we measured the 
12 C(p, n) reaction at 9\ a b = 0° and incident beam energies of T p = 198 and 
295 MeV. The target thickness is 38 mg/cm 2 whose contribution to the final 
energy resolution is negligible small. The results are shown in Fig. 5. The 
energy resolutions were evaluated by fitting the peak of the 12 N state to the 
standard hyper-Gaussian and they are about 280 and 500 keV FWHM for T p 
= 198 and 295 MeV data, respectively. These energy resolutions are consistent 
with the simulation results shown in Fig. 3, and they satisfy our requirements. 
These counters have been used in the NPOL3 which will discussed in the next 
section. 



3 Neutron detector NPOL3 and its performance 

3. 1 Neutron detector NPOL3 

The neutron detector NPOL3 has two main neutron detector planes, HD1 
and HD2. Each plane consists of 10 stacked 100 cm x 10 cm x 5cm plastic 
scintillator bars made of Bicron BC408. Two thin plastic scintillator planes 
CPV and CPD, used to identify charged particles, complete the NPOL3 which 
is schematically shown in Fig. 4. Additional two-dimensional neutron counter 
NC is used in polarimetry mode, which will be discussed in Sec. 4. A scintillator 
in HD1 or HD2 is viewed by the fast Hamamatsu 25] H2431 (HD1) or H1949 
(HD2) PMTs attached at each end. The one-dimensional position is deduced 
from the fast timing information derived from PMT. The CPV and CPD are 
made of three sets of plastic scintillators (BC408) with a size of 102 cm x 35 
cm x 0.5 cm. The scintillator is viewed from both sides (U:up and D:down) by 
the Hamamatsu H7195 PMT through fish-tail shape light guides. The CPV is 
used to veto charged particles entered into the NPOL3, and the CPD is used 
to detect the recoiled charged particles. 

In the neutron detector mode of NPOL3, a neutron is detected by either 
neutron detector plane HD1 or HD2. Furthermore we have required that the 
recoiled charged particle should be detected by the following detector plane. 
This means that we measure neutrons in the polarimetry (double-scattering) 
mode of the NPOL3 system. This procedure has been applied in the data 
analysis of the INPOL system j^], and it is useful to improve the final energy 
resolution as well as to reduce the low-energy tail component. Thus we have 
applied this procedure in our data analysis. 
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3.2 Light output calibration 



The neutron detection efficiency depends on the threshold level £ t h of the 
scintillator light output. The total light output £ is constructed by taking the 
geometric mean of the PMT outputs at both ends, Cl and jCr, as 



The light output C is calibrated by using 4.4 MeV 7-rays from a 241 Am- 9 Be 
neutron source and cosmic rays (mostly fi^). In the cosmic-ray measurement, 
10-bar hits in either HD1 or HD2 were required. Thus horizontal and zenith 
angles can be measured in HD1 and HD2, respectively, on the assumption 
that cosmic rays pass in straight lines through the detector plane. The £ for 
cosmic rays has been corrected for these angles in order to evaluate the C in 
case cosmic rays pass the scintillator of 10 cm in thickness. 

The left and right panels of Fig. 6 show the C spectra for 4.4 MeV 7-rays and 
cosmic rays, respectively. We can clearly see the Compton edge (4.2 MeV cc ) for 
the 7-ray spectrum and the sharp peak with a Landau tail for the cosmic-ray 
spectrum. The histograms represent the results of the Monte-Carlo simula- 
tions with geant4 described in Sec. 2. The simulations successfully repro- 
duce the main components of the Compton edge and the sharp peak in 7-ray 
and cosmic-ray spectra, respectively. Furthermore they well describe the one- 
photon escape peak at ~3.5 MeV ee in the 7-ray spectrum and the Landau 
tail in the cosmic-ray spectrum. Note that the discrepancy between the data 
and the simulation result for 4.4 MeV 7-rays in the low C region is due to the 
contributions from 7-rays other than 4.4 MeV 7-rays from the neutron source 
which are not taken into account in the simulations. 

3.3 Neutron detection efficiency 

The laboratory differential cross section da^/dVL is related to the number of 
measured neutrons N„ as 



where N p is the number of incident protons, p is the target density, AQ is the 
solid angle, e is the intrinsic neutron detection efficiency, T is the transmission 
factor along the neutron flight path, and fu ve is the experimental live time. 

We measured the the product eT by measuring the neutron yield from the 0° 
7 Li(p, n) 7 Be(g.s. + 0.43 MeV) reaction which has a constant center of mass 




(4) 



daiab N n 



(5) 



dVL " N pP AneTf hvc ' 
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(cm.) cross section of a r rn (0°) = 27.0 ± 0.8 mb/sr at the bombarding energy 
range of T p = 80-795 MeV j27J. The measurements were performed at T p = 198 
and 295 MeV by using the enriched 7 Li target with a thickness of 54 mg/cm 2 . 
The integrated beam current N p was measured by using an electrically isolated 
graphite beam stop (Faraday cup) connected to a current integrator. The live 
time /ii ve was typically 90%. 

The finally obtained neutron detection efficiencies (eT) are shown in Fig. 7 
as a function of the light output threshold Ah- The systematic uncertainty 
is estimated to be about 6% by taking into account of the uncertainties both 
of the 7 Li(p, n) cross section and of the 7 Li target thickness. It is found that 
the detection efficiencies are almost independent of the neutron kinetic energy 
with a value of ~0.017 at £ t h = 5 MeV ee - 



4 Neutron polarimeter NPOL3 and calibration procedure 

4-1 Neutron polarimeter NPOL3 

In the polarimetry mode of NPOL3, the neutron detector planes, HD1 and 
HD2, serve as neutron polarization analyzers, and the following two-dimensional 
position-sensitive neutron counter NC acts as a catcher of doubly scattered 
neutrons or recoiled protons. Neutron polarization is determined from the 
asymmetry of the n + p events whose analyzing powers can be rigorously ob- 
tained by using the phase-shift analysis of the NN scattering. The n+p events 
are kinematically resolved from the other events by using time, position, and 
pulse-height information from both analyzer and catcher planes. Note that 
this kinematical selection significantly reduces the background events from 
the wraparound of slow neutrons from preceding beam pulses, cosmic rays, 
and the target gamma rays. However, the quasielastic n + C reaction has an 
influence on the determination of the neutron polarization because its kine- 
matical condition is similar to that of the n + p scattering. Thus we have 
measured the effective analyzing powers which include the contributions from 
both n + p and n + C by using the polarized neutron beams. 

4 ■ 2 Effective analyzing powers 

There are two channels, (n,n) and (n,p), in the polarimetry mode. Doubly 
scattered neutrons and recoiled protons are measured by the catcher plane 
NC in (n, n) and (n, p) channels, respectively. The effective analyzing powers, 
A™ s and Ay^ eS , of (n, n) and (n, p) channels can be determined by using the 
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neutron beam with a known polarization p' N . Ay™ s and Ay? eS are deduced by 

^nn ^np 

A y?es = A y;eS = ~r> ( 6 ) 

Pn Pn 

where e nn and e np are the asymmetries measured for (n, n) and (n, p) channels, 
respectively. 



4-3 Polarized neutron beams 



The first source of polarized neutrons is the 0° 2 H(p, n)pp reaction at T p = 
198 MeV. The Dll value at T p = 198 MeV was obtained by using the data 
in Ref. j^. Because this reaction is mainly a Gamow- Teller transition, the 
polarization transfer coefficients D{j satisfy ^| 

The neutron polarization p' N is given by using Dnn(0°) as 

p' N = p N D NN (0°), (8) 
where p^ is the polarization of the incident proton beam. 

The second source is the 0° 6 Li(p, ra) 6 Be(g.s.) reaction at T p = 198 MeV. This 
reaction is also a GT transition, and its D NN value was measured at T p = 
200 MeV by Taddeucci 29]. The result can be used as a double-check of the 
calibration performed by using the neutrons from 2 H(p, n)pp. 

The third source is the 0° 12 C(p, n) 12 N(g.s.) reaction at T p = 198 MeV. The 
Dnn{0°) value of this GT transition was also measured at T p = 200 MeV by 
Taddeucci j2^|. Unfortunately, its uncertainty is relatively large compared with 
that of the 2 H(p, n) or 6 Li(jo, n) 6 Be(g.s.) reaction. Thus the final accuracy of 
Ay- e s is limited by the uncertainty of the neutron beam polarization originating 
from the uncertainty of Dnn{0°), Therefore we have used the following method 
which does not require to know the Da(0°) values. 

The GT transition is the spin-flip and unnatural-parity transition, and its po- 
larization transfer coefficients satisfy the relation in Eq. (7). Here we assume 
that two kinds of polarized proton beams are available; one has pure longi- 
tudinal polarization pi and the other has pure normal polarization p^- The 
neutron polarizations at 0° become p' L = PlDll(0°) and p' N = pnDnn(0°) for 
the beam polarizations of pl and p^, respectively. Then the resulting asym- 
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metries measured by a neutron polarimeter are 



e L = p' L A y;cS = p L D LL (0 o )A y . cS , 

e N = PN^-y^S = PND NN (0°)A yie s. 



(9) 



By using Eqs. (7) and (9), Ay- e s can be expressed as 



\Pl PnJ 



(10) 



Thus the A y . c g value can be obtained without knowing a priori the values of 
Du(0°), and we have applied this technique to the 12 C(p, n) 12 N(g.s.) reaction. 

4-4 Experimental conditions 

We used a deuterated polyethylene (CD 2 ) target with a thickness of 228 
mg/cm 2 as a deuteron target. Neutrons from the 2 H(p, n)pp reaction were 
clearly separated from those from the 12 C(p, nX) reaction since the former 
reaction Q value of —2.2 MeV is significantly smaller than the latter Q value 
of > —18.1 MeV. We also used 99% enriched 6 Li and natural C targets with 
thicknesses of 181 and 38 mg/cm 2 , respectively. The proton beam intensity 
and its polarization are typically 80 nA and 0.70, respectively. 

4-5 Sector methods 

We have adopted the sector method in order to obtain the asymmetry. In 
this method, the plane (9, <f>) defined by the scattering polar angle 9 and the 
azimuthal angle is divided into four sectors: Left, Right, Up, and Down, as 
is shown in Fig. 8. Each sector is defined by scattering angle limits 9 mm and 
Qmax an( j j- ne az j mu thal half-angle $. 

The asymmetry e is calculated by using the event numbers in these sectors. 
For example, the relation between the left-right asymmetry e'^ and the event 
numbers is as follows. The azimuthal distribution of scattered particles can be 
described as 



where we suppress the 9 dependence for simplicity. Then the event numbers 
scattered to left and right become 



(11) 




N 



7T 




(12) 




N 




7T 
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where N is the event number accepted in the 9 bin, and the acceptance is 
assumed to be constant. Thus the asymmetry = p' N A y . ef[ is found to be 



e> - Yl ~ Yr $ (13) 

% -y L + ^sin$- {l6) 

Note that the sideway polarization p' s is deduced from the up-down asymmetry 
e' s which is obtained by substituting L — > D and R — ► U in Eq. (13). 



5 Optimization of performance of NPOL3 



5.1 Optimization criteria 



The performance of a polarimeter can be measured by its figure-of-merit 
(FOM) value defined as 

FOM = N(A y , c{i ) 2 , (14) 

where N is the total event number Yl + Yr or Yjj + Yd. We can minimize the 
uncertainty of the final result (neutron polarization and Dij) by maximizing 
this FOM value. 

We have applied three types (Types-I, II, and III) of software cuts to the mea- 
sured 9, 0, inter-plane velocity, and pulse height information of the NPOL3. 
The FOM value is maximized in Type-I, whereas the A y . cS is maximized in 
Type-II with keeping the FOM value larger than 40% of the optimum value. 
In Type-Ill, both FOM and A y . e s values have been optimized with keeping 
the FOM value larger than 80% of the optimum value. The Type-I is the best 
choice from the statistical point of view in order to maximize the performance 
of the NPOL3. However, in general, a higher A y . eS value in Type-II or III is 
helpful to reduce the systematic uncertainty of the final result. Thus we have 
obtained these three sets and we will select the optimum set depending on the 
experimental condition. 

In the following, we show the results by using the data of the 2 H(j9, n)pp 
reaction because the total event number is significantly larger than that of 
other two polarized-neutron production reactions. The A y . c g value is also ob- 
tained from the 6 Li(p, n) 6 Be(g.s.) data by applying the same software cuts. 
The mean neutron energy is about 193 MeV for both the the 2 H(p, n)pp 
and 6 Li(p, n) 6 Be(g.s.) reactions. Thus we can use the results for a double 
check of the calibration. The same software cuts are also applied to the 
12 C(j9, n) 12 N(g.s.) data. The mean neutron energy of ~180 MeV is signifi- 
cantly low, therefore, the result can be used to discuss the neutron energy 
dependence of A y - C s. 
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5.2 Optimization o/$ 



Figure 9 shows the A y . e g and FOM values in Type-I as a function of the 
azimuthal half- angle $. The FOM values are normalized to a maximum FOM 
value of 1. The A y;e s values are insensitive to the choice of $, whereas the 
FOM values take a maximum at $ ~ 67°. This angle is consistent with the 



simple estimation of the optimum $ = 66.8° |2jj. Because the $ dependence 



is almost independent of the software cut types, we set $ = 66.8° in all cases. 



5.3 Optimization of 9 window 



The event selection on 9 is necessary to prevent dilution of A y . e g by events 
with small or negative analyzing powers. Figure 10 shows the scattered-particle 
distribution as a function of 9 in Type-I. The results in other Types-II and 
III are also identical. The double scattering events with a scattering angle up 
to ~ 60° are accepted by the NPOL3. The analyzing power A y {9) of the free 
n+p scattering is positive only up to 45° and 41° for (n, n) and (n,p) channels, 
respectively, at T n = 193 MeV. Thus part of the events should be rejected on 
the basis of small A y as described below. 

Figure 11 shows the contour plots of the FOM values in Type-I as functions 
of both 9 mm and # max . Figure 12 represents the A y . e s and FOM values as a 
function of either 9 mm or # max . In both figures, the FOM values are normalized 
to a maximum FOM value of 1. The FOM values strongly depend on both 9 mm 
and # max ; whereas the A y . e g values weakly depend on them except for # max in 
the (n, n) channel. In cases of Types-II and III, the optimal (9 mm , # max ) values 
are selected on basis of the criteria described in Sec. 5.1. In Table 1, we present 
the finally chosen (9 min , # max ) values. 



5.4 Optimization of R v window 



The kinematical selection of the n + p events is performed by using the kine- 
matical quantity R v defined as 

Rv = — , (15) 

VNN 

where v is the measured particle velocity between the analyzer and catcher 
planes and vjyjy is the predicted velocity based on the NN kinematics. The 
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vnn value is given by 



2cos6JT n /(T n + 2m N ) 

V NN = C- r, 16 

1 - cos 9 T n j [T n + m N ) 

where c is the light velocity and is the nucleon mass. Figure 13 shows 
the scattered-particle distribution in Type-I as a function of R v . The peaks at 
R v ~ 1 correspond to the events from the p(n,n)p or p(n,p)n reaction. The 
shoulders at < 1 are due to the contributions from quasielastic scattering, 
inelastic scattering, and wraparounds. 

Figure 14 shows the contour plots of the FOM values as functions of both 
i?™ in and i?™ ax . Figure 15 represents the A y;e g and FOM values as a function 
of either i2™ m or i?™ ax . The FOM values are normalized to a maximum FOM 
value of 1 in both figures. It is found that the FOM values are almost insensitive 
to i?™ ax for i?™ ax > 1.3. However, i?™ ax has been limited to be less than 1.3 
to eliminated the contribution from 7-ray events with R v > 1.5. In Types-II 
and III, the optimal (-R™ m , -R™ ax ) values are selected on basis of the criteria 
described in Sec. 5.1. The final (-R™ in , -R™ ax ) values chosen are listed in Table 1. 



5.5 Optimization o/£ th 



Figure 16 shows the A y . eS and FOM values in Type-I as a function of C t h- 
The FOM values monotonously decrease as increasing C th , whereas the A y . e s 
values are almost independent of £ t h- In all software cut types, we have chosen 
£ t h = 4 MeV ee to ensure reproducibility in the offline analysis rather than just 
accept the hardware threshold of ~2 MeV ee . 



6 Performance of polarimeter NPOL3 



6. 1 Effective analyzing powers 



The effective analyzing powers A y . ea obtained from the 2 H(p, n)pp data by 
applying previously described software cuts are summarized in Table 2. The 
A y]e s values deduced from the 6 Li(p, n) 6 Be(g.s.) data are also listed in Table 2. 
The results of the 2 H(p, n)pp and 6 Li(p, n) 6 Be(g.s.) data are consistent with 
each other within their statistical uncertainties, which suggests the reliability 
of the present calibration procedure. 

The A y . ef [ values deduced from the 12 C(p, n) 12 N(g.s.) data are plotted in 
Fig. 17, and they are also listed in Table 2. The results of the 2 H(p, n)pp 
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data are also displayed in Fig. 17. The A y - e g values have been estimated as a 
function of T n on the basis of the free n + p analyzing powers derived from 
the NN phase shift analysis j^. In the calculations, the (9,4>) cuts and the 
counter geometries have been properly taken into account. The solid curves 
in Fig. 17 show the results of the calculations. The measured A y;e s values are 
different from the calculated values by factors of 0.73-0.79 and 1.16-1.29 for 
(ft, n) and (n,p) channels, respectively, depending on the event cut types. The 
difference of A y . e g is mainly due to the contributions from the quasielastic 
n + C events. The dashed curves represent the calculated A y . eS values normal- 
ized to reproduced the experimental results. These normalized values will be 
used in future data analysis in case there is no appropriate calibration result. 



6.2 Double scattering efficiency 

The double scattering efficiencies (cd.s.) of the NPOL3 were deduced from 
the data of the 2 H(p, n)pp and 12 C(p, n) 12 N(g.s.) reactions. The results are 
listed in Table 3. The 6d.s. values in both (ft, n) and (n,p) channels are almost 
constant in the present neutron energy region. 



7 Discussion 



7.1 9 dependence of A y . eS 

In the sector method, the events within (9 mm , # max ) are integrated in order to 
evaluate the A y . C Q values in the calibration. Thus the results can be considered 
as the mean effective analyzing powers in the range of 9 = # mm -# max . We 
deduced the A y . eS values as a function of 9 in order to check whether the 
n + p events are properly selected in our analysis. As a typical example, the 
results in Type-Ill are shown in Fig. 18. The solid curves display the angular 
distributions of A y (9) of the free n + p scattering. The measured angular 
distributions are very similar to those of the free n+p scattering. However the 
magnitudes are different mainly due to the contributions from the quasielastic 
n + C events. The dashed curves are the results normalized to reproduce 
the measured values, and the normalization factors are 0.7 and 1.2 for (n, n) 
and (n, p) channels, respectively. The agreement of the angular distributions 
supports the proper kinematical selection of the n + p events in our analysis. 
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7.2 Comparison with INPOL system 



Here we compare the performance of the NPOL3 system with that of the 
INPOL system [26| which has been used in the same energy region. For the 
(n, n) channel, the NPOL3 has ~l/2 volume for both analyzer and catcher 
planes compared with the INPOL. Thus the double scattering efficiency cd.s 
and the corresponding FOM value of the NPOL3 is expected to be ~l/4 of 
those of the INPOL. The FOM value of the NPOL3 is 0.28 x 1(T 4 in Type- 
I at T n = 193 MeV, which is about 1/4 of that of the INPOL reported as 
0.92 x 10~ 4 at T n = 194 MeV. For the (n,p) channel, the FOM value of the 
NPOL3 is expected to be similar to that of the INPOL because their volume 
of the analyzer planes is similar with each other. In fact, the FOM value of 
0.89 x 10~ 4 in Type-I is close to that of 0.81 x 10~ 4 of the INPOL. Thus we 
conclude that the performance of the NPOL3 is very similar to that of the 
INPOL by taking into account of the difference of the thicknesses of both 
analyzer and catcher planes. 



8 Summary 



The high resolution neutron polarimeter NPOL3 has been constructed and 
developed to measure polarization transfer observables Dij for (p, n) reactions 
at intermediate energies around T p = 200 MeV. The NPOL3 system can mea- 
sure the normal as well as sideways components of the neutron polarization 
simultaneously. The other longitudinal component can be measured by using 
the neutron spin rotation magnet described in Ref. ji^. Thus we can perform 
the measurement of a complete set of with a high resolution of ~300 keV 
by using the NPOL3. 

The performance of the NPOL3 system was studied by using the polarized 
neutrons from 2 H(p, n)pp, 6 Li(p, n) 6 Be(g.s.), and l2 C(p, n) 12 N(g.s.) reactions 
at T p = 198 MeV. The effective analyzing powers in Type-I are 0.33 and 0.12 
for (ft, n) and (n,p) channels, respectively. The performance is comparable to 
that of the INPOL system at the same energy. 
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Fig. 1. Illustration of the Monte-Carlo simulation with geant4. The neutron beam 
with T n = 200 or 300 MeV bombard a plastic scintillator with a size of 100 cm x 10 
cm x AL cm. The optical photons generated by recoiled charge particles propagate 
in the scintillator and are detected by PMTs at both ends. The time information, 
ti and tR, deduced from PMT signals are used to evaluate the time of flight t and 
its width At. 

Fig. 2. The pulse-height (upper panel) and CFD timing (lower panel) spectra 
generated by the Monte-Carlo simulation with geant4. See text for details. 

Fig. 3. The simulation results for the energy resolution AT n as a function of the 
counter thickness AL. The filled circles and the filled boxes represent the results for 
the measurement at T p = 200 and 300 MeV, respectively, as explained in the text. 

Fig. 4. Schematic (left) and top (right) views of the NPOL3 system. In the po- 
larimetry mode of NPOL3, HD1 and HD2 are the analyzer planes while NC is the 
catcher plane. Thin plastic scintillator planes are used to veto (CPV) or identify 
(CPD) charged particles. 

Fig. 5. Measurements of the 12 C(p, n) reaction at T p = 198 MeV (upper panel) and 
297 (lower panel) MeV with 5 cm-thick plastic scintillation counters. The curve of 
each panel shows the reproduction of the 12 N(g.s.) state with the hyper-Gaussian 
peak. 

Fig. 6. Light output spectra for 4.4 MeV 7-rays from a 241 Am- 9 Be neutron source 
(left panel) and cosmic rays (right panel). The histograms represent the results of 
the Monte-Carlo simulations with geant4. See text for details. 

Fig. 7. Neutron detection efficiencies of 5-cm thick plastic scintillation counters as a 
function of detection threshold £ t h- The filled circles and the filled boxes represent 
the results at T n = 190 and 291 MeV, respectively. 

Fig. 8. Illustration of the sector definition on the catcher plane NC for left/right 
(left panel) and up/down (right panel) asymmetries. 

Fig. 9. Effective analyzing powers A y . e Q and FOM as a function of azimuthal 
half-angle <!>. The upper and lower panels correspond to the results for (n, n) and 
(n,p) channels, respectively. 

Fig. 10. Distributions of the doubly scattered particles as a function of scattering 
polar angle 6. The left and right panels correspond to the results for (n, n) and 
(n,p) channels, respectively. 

Fig. 11. FOM of NPOL3 as functions of lower and upper limits of the scattering 
polar angle window (g mm ^0 max y The left and right panels correspond to the re- 
sults for (n, n) and (ft, p) channels, respectively. The selected limits in Type-I are 
indicated by the filled circles. 
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Fig. 12. Effective analyzing powers A y;e fi and FOM as a function of either lower (left 
panels) or upper (right panels) limit of the scattering polar angle 6. The upper and 
lower panels correspond to the results for (n,n) and (n,p) channels, respectively. 



Fig. 13. Distributions of the doubly scattered particles as a function of velocity 
ratio R v . The left and right panels correspond to the results for (n, n) and (n,p) 
channels, respectively. 



Fig. 14. FOM of NPOL3 as functions of lower and upper limits of the velocity ratio 
window (-R™ in , i?™ ax ). The left and right panels correspond to the results for (ft, n) 
and (n,p) channels, respectively. The selected limits in Type-I are indicated by the 
filled circles. 



Fig. 15. Effective analyzing powers A y . e Q and FOM as a function of either lower 
(left panels) or upper (right panels) limit of the velocity ratio R v . The upper and 
lower panels correspond to the results for (n, n) and (n,p) channels, respectively. 



Fig. 16. Effective analyzing powers A y - C g and FOM as a function of detection 
threshold C t h- The upper and lower panels correspond to the results for (n, n) and 
(n,p) channels, respectively. 



Fig. 17. Effective analyzing powers of the NPOL3 system for (n, n) (filled circles) 
and (n,p) (filled boxes) channels, respectively, as a function of neutron energy. The 
top, middle, and bottom panels correspond to the results in the software cut Types-I, 
II, and III, respectively. The solid curves represent the calculations based on the 
free n+p observables. The dashed curves represent the calculations normalized to 
reproduce the experimental data. 



Fig. 18. Angular distributions of effective analyzing powers of the NPOL3 system 
for (n, n) (filled circles) and (n,p) (filled boxes) channels, respectively. The solid 
curves represent the calculations based on the free n + p observables whereas the 
dashed curve represent the normalized results to reproduce the experimental data. 
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Type-I Type-II Type-Ill 

Quantity FOM optimum A y . e ff optimum Intermediate 





(n, n) 


(n,p) 


(n, n) 


(n,p) 


(re, n) 


(re» 




67° 


67° 


67° 


67° 


67° 


67° 


8 


(8°, 30°) 


(8°, 53°) 


(12°, 24°) 


(20°, 53°) 


(10°, 27°) 


(12°, 53°) 


R v 


(0.90,1.30) 


(0.70, 1.30) 


(0.96,1.30) 


(0.88,1.30) 


(0.92,1.30) 


(0.70,1.30) 


Ah 


4MeV ee 


4MeV ee 


4MeV ee 


4MeV ee 


4MeV ee 


4MeV ee 



Table 1 



Optimized cuts used in the NPOL3 system. 



Type-I Type-II Type-Ill 

Energy (Reaction) FOM optimum A y . e fi optimum Intermediate 





(re, n) 


(n,p) 


(n, n) 


(n,p) 


(re, n) 


(n,p) 


193 MeV ( 2 H(p, n)pp) 


0.333 


0.123 


0.398 


0.154 


0.373 


0.128 




±0.008 


±0.002 


±0.013 


±0.003 


±0.009 


±0.002 


193 MeV ( 6 Li(p,n) 6 Be(g.s.)) 


0.305 


0.123 


0.369 


0.148 


0.341 


0.136 




±0.025 


±0.009 


±0.032 


±0.011 


±0.028 


±0.010 


180 MeV ( 12 C(p,ra) 12 N(g.s.)) 


0.358 


0.105 


0.427 


0.129 


0.398 


0.112 




±0.021 


±0.004 


±0.034 


±0.008 


±0.024 


±0.005 



Table 2 



Effective analyzing powers of the NPOL3 system. 



Type-I Type-II Type-Ill 

Energy FOM optimum ^4y ;c ff optimum Intermediate 

(n,n) (n,p) (n,n) (n,p) (ft,n) (n,p) 

193 MeV 2.5 x 10~ 4 5.9 x 10~ 3 0.9 x 10~ 4 1.7 x 10~ 3 1.8 x 10~ 4 4.9 x 10~ 3 

180 MeV 2.8 x 10~ 4 6.5 x 10~ 3 1.0 x 10~ 4 1.8 x 10" 3 2.0 x 10~ 4 5.4 x 10" 3 
Table 3 

Double scattering efficiencies of the NPOL3 system. 
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